Relativistic Outflow of Electron-Positron Pair Plasma from a 

Wien Equilibrium State 
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ABSTRACT 

We investigate a novel mechanism of bulk acceleration of relativistic outflows of pure electron 
positron pairs both analytically and numerically. The steady and spherically symmetric flow is 
assumed to start from a Wien equilibrium state between pure pairs and photons in a compact 
region which is optically thick to the electron scattering at a relativistic temperature. Inside the 
photosphere where the optical thickness becomes unity, pairs and photons behave as a single fluid 
and are thermally accelerated. Outside the photosphere, pairs and photons behave separately 
and we assume the free streaming approximation for photons, which are emitted from the rela- 
tivistically moving photosphere. Pairs are shown to be thermally accelerated further even outside 
the photosphere because the photospheric temperature is at least mildly relativistic. It is to be 
noted that the mean energy of photons is higher than that of pairs in the comoving frame of pairs 
and that the Compton interaction leads to additional heating and radiative acceleration of pairs. 
For a reasonable range of the boundary temperature and optical thickness, the terminal Lorentz 
factor of pair outflows turns out be more than 10 and the terminal kinetic power accounts for 
more than 2/3 of the total luminosity. While the total luminosity should be at least larger than 
the Eddington luminosity, real luminosity can be modest if the outflow is collimated by some 
unknown mechanisms. This mechanism successfully avoids the difficulties of pair annihilation 
and radiation drag owing to the pair production by accompanying high energy photons and the 
strong beaming of radiation field. It is seen that most pairs injected at the boundary survive 
to infinity. The radiation from the photosphere should be observed as MeV peaked emission at 
infinity with an order of kinetic power of jets. 

Subject headings: electron-positron pairs, hydrodynamics, relativistic plasma, relativistic jets, blazars 



1. Introduction 

The production and bulk acceleration of rela- 
tivistic jets in active galactic nuclei is one of the 
most challenging problems in astrophysics. The 
basic properties to be explained are the highly 
relativistic velocity with the bulk Lorentz factor 
of around 10, the kinetic power almost compara- 
ble to the Eddington luminosity, and the collima- 
tion into a small opening angle (e.g., Ostrowski 
et al. 1997; Begelman, Blandford and Rees 1984). 
Although many ideas have been proposed rang- 



ing from magnetic and radiative accelerations to 
Blandford-Znajek process, each model has its own 
difficulties and there is no general consensus on 
how jets are produced and accelerated (e.g., Begel- 
man 1995; Begelman, Blandford and Rees 1984). 
In this paper, we explore a thermal acceleration 
mechanism of pure electron-positron plasma in 
view of the recent observational and theoretical 
developments. 

Several important observational developments 
have been made in this decade. First of all, 
EGRET on board the Compton Gamma Ray Ob- 
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servatory has found that blazars are a strong 7-ray 
emitter (Mukherjee et al. 1997). The 7-ray emis- 
sion is interpreted in terms of the inverse Compton 
scattering off relativistic electrons in the jets the 
direction of which is so close to the line of sight 
that the emission from jets is strongly beamed 
(e.g., Blandford & Konigl 1979; Maraschi, Ghis- 
cllini and Celotti 1992; Sikora, Begelman and Rees 
1994). Combining the synchrotron emission at 
lower energies, we can estimate the physical states 
of the emission region in the jets. The results of 
multi-frequency spectral fitting indicate that the 
jets are particle dominated rather than magnetic 
field dominated (Inoue and Takahara 1996; Kino, 
Takahara and Kusunose 2001). The energy den- 
sity of relativistic electrons is more than one or- 
der of magnitude larger than that of the magnetic 
field. The size of the emission region is estimated 
to be 10 16 ~ 10 17 cm from time variabilities and 
visibility of 7-ray photons. Relativistic beaming 
factor of around 10 is required to assure the vis- 
ibility of 7-ray photons. Since the distance of 
the emission region from the central black hole 
is reasonably estimated to be 10 times its size, it 
is around 10 17 ~ 10 18 cm, an order of magnitude 
shorter than that of supcrluminal knots observed 
at radio wavelengths. Thus, the bulk acceleration 
of jets must occur within such short distances. As 
for the matter content of jets, electron-positron 
pair dominance has been proposed based on sev- 
eral independent observations (Wardle et al. 1998; 
Reynolds et al. 1996; Hirotani et al. 1999, 2001), 
although the contribution of protons is not com- 
pletely excluded. 

Second important development is the estima- 
tion of the kinetic power of extended radio sources. 
Such studies have shown that large scale kinetic 
power is well correlated with radiative power and 
suggest that jets are produced by accretion process 
rather than they are two unrelated phenomena 
(Rawlings and Saunders 1991). Related to this, 
Galactic super luminal sources GRO J 1655+40 and 
GRS 1915+105 with the bulk Lorentz factor of 
a few were discovered (Hjcllming & Rupen 1995; 
Mirabel & Rodriguez 1994). The jet ejection for 
these two sources seems to be correlated with X- 
ray behavior, which suggests that the jets are re- 
ally produced by accretion (Feroci et al. 1999). 

Finally, developments of observations and mod- 
els of cosmic gamma-ray bursts (GRBs) are quite 



fascinating. Confirmation of the afterglows of 
GRBs implies the basic correctness of the fireball 
model (Piran, Shemi and Narayan 1993; Rees and 
Meszaros 1992). In this model, GRBs are ini- 
tiated by the production of a compact high en- 
tropy fireball, although no specific mechanism of 
the production process itself has been widely ac- 
cepted. The fireball is relativistic thermal plasma, 
which consists of photons, electron-positron pairs 
and a small amount of baryons. As the fireball ex- 
pands, it is thermally accelerated and the temper- 
ature drops. When the temperature drops below 
the electron mass, electron-positron pairs annihi- 
late, but it is still optically thick to the electron 
scattering because of a small baryon load. Final 
outcome is that the initial thermal energy is con- 
verted to the bulk kinetic energy of baryons at- 
taining the bulk Lorentz factor of E/M^c 2 , where 
E and Mb are the total energy and baryon mass 
of the fireball. Inhomogeneitics in the fireball are 
expected to produce internal shocks which are sup- 
posed to correspond to GRBs themselves and an 
impact to ambient matter produces an external 
shock which corresponds to afterglows. The emis- 
sion from shocks is ascribed to the particle accel- 
eration by shocks. It is also mentioned that some 
GRBs are described by jets rather than spherical 
explosion with respect to the afterglow light curve 
and energetics. 

We believe that the success of the fireball 
model for GRBs has strong implications for jets 
in blazars. For blazars, internal shocks corre- 
spond to the compact emission regions which are 
located relatively near the central black hole and 
external shocks produce large-scale radio sources 
(Carilli and Barthel 1996; Begelman and Cioffi 
1989). But, one strong difference between blazars 
and GRBs is the size of the fireball. The initial 
size of GRB fireballs is 10 6 ~ 10 7 cm and it is in 
a complete thermal equilibrium for the total en- 
ergy of ~ 10 52 erg. In contrast, the size of blazar 
fireballs should be around 3r g <~ 10 14 cm with a 
similar total energy, where r g = 2GM/c 2 is the 
Schwarzshild radius for the black hole mass M; 
the typical value of M is taken to be 10 8 M©. If 
the blazar fireball is in a complete equilibrium, 
the temperature is an order of 10 5 K which is far 
below the electron mass and no pair production is 
expected. Basically it describes an optical thick 
radiation pressure dominated standard accretion 



2 



disk and no relativistic expansion is expected. 

However, for an optically thin state, elec- 
tron temperature can be relativistic and copi- 
ous electron-positron pairs can exist in principle. 
Even for mildly relativistic temperatures, copious 
electron-positron pairs exist for a finite optical 
thickness to scattering (Svensson 1984; Lightman 
1982; Kusunose and Takahara 1985). It has been 
shown that the pair density is limited to a rela- 
tively small value in hot accretion disks when the 
pair production is balanced with pair annihilation 
(Kusunose and Takahara 1988). However, effect of 
pair escape can increase the pair density as was ex- 
emplified by Yamasaki, Takahara, and Kusunose 
(1999). They have argued that pair outflow from 
hot accretion disks may account for major part of 
the dissipated energy from accretion. In such a 
state, photons and electron-positron pairs are not 
in a complete thermal equilibrium but they are 
coupled by electron scattering and pair produc- 
tion/annihilation process and the photon density 
is far below that of the equilibrium state for a 
given temperature. Thus, the escaping pairs just 
above the disk may act as the fireball to produce 
relativistic jets. Although the physical state of 
such escaping pairs has not been fully investi- 
gated by now, we believe that this is the most 
promising way of producing relativistic jets with 
electron-positron pairs. 

In this paper, we consider the bulk acceleration 
of such fireballs starting from a range of the initial 
states in the simplest way. We treat steady and 
spherically symmetric outflows starting from pure 
electron-positron pairs with photons in a Wien 
equilibrium for which photon distribution function 
is given by a Boltzmann distribution with a finite 
chemical potential. The formulation is similar to 
that of Grimsrud & Wasserman (1998) (GW98), 
in which relativistic flows of pure electron-positron 
pairs starting from a complete thermal equilib- 
rium were analyzed. The dynamics starting from 
a Wien equilibrium is expected to be similar to 
that given by GW98 when proper accounts for 
the difference arc taken. This type of flows takes 
into account of the effects of pair production as 
well as pair annihilation. Then, we can examine 
to what extent we can avoid the pair annihila- 
tion problem pointed out by some authors arguing 
against electron-positron pair dominance in rela- 
tivistic jets (Celotti and Fabian 1993; Blandford 



and Levinson 1995). If pairs eventually annihilate 
in the course of expansion, there would remain 
very little kinetic power. It is very important to 
quantitatively assess the surviving kinetic power 
of such a fireball. 

In section 2, we give basic formulation, in sec- 
tion 3 we present approximate analytical solutions 
and in section 4 we present the results of numerical 
solutions. 

2. Formulation 

We assume that at the inner boundary the 
source of pure electron-positron pairs and high en- 
ergy photons injects a steady outflow of pairs and 
radiation. The optical thickness of such a region 
to the electron scattering is large but that to the 
absorption is smaller than unity. Thus, the phys- 
ical state of such a regime may be characterized 
by a Wien equilibrium in the simplest case. When 
we neglect the dynamical effects, physical state of 
such a plasma with a finite amount of protons is 
known to be strongly constrained by thermal bal- 
ance and pair equilibrium owing to the inefficiency 
of pair annihilation at relativistic temperatures 
(Lightman 1982; Svensson 1984; Kusunose and 
Takahara 1985). However, we should note that 
the time scale to reach such a steady state is fairly 
long. If the plasma is first put at a relativistic 
temperature without pairs, pair production pro- 
ceeds very quickly and the subsequent photon pro- 
duction and Compton scattering strongly cool the 
plasma as the pair concentration still increases. 
Finally, pair annihilation leads to a pair equilib- 
rium state. The time scale to reach the pair equi- 
librium state is controlled by the photon produc- 
tion and photon escape, the time scales of which 
are much longer than the dynamical time scale. 
Thus, the dynamical effects are essential to the 
behavior of such plasma and most probably lead 
to the expansion and outflow. The constraints 
on the electron temperature and optical thickness 
should be drastically changed from those for static 
confined plasma. In this paper, for simplicity, we 
assume that a Wien equilibrium state is already 
established from the outset and that the tempera- 
ture and optical thickness can take a wider range 
than those of static cases. Since the properties 
of the initial states, where dynamical effects are 
critically important, are not well known, we here 
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regard the temperature and optical thickness at 
the boundary as free parameters. 

The flow is optically thick near the bound- 
ary and becomes optically thin beyond the pho- 
tosphere. We treat separately optically thick and 
thin regimes as GW98. For optically thick regime 
in which the optical thickness to the electron 
scattering for radially outgoing photons is larger 
than unity, we treat as if electron-positron pairs 
and photons form a single fluid. Although, ra- 
diation field may be anisotropic, we simply as- 
sume isotropic photon distribution in the comov- 
ing frame. Note that as was explicitly shown in 
GW98, the isotropy of photon field is a surpris- 
ingly good assumption in a relativistic expanding 
flow. Although the assumption of a Wicn equi- 
librium is of course a simplification to treat the 
dynamical problem in a simple way, we may re- 
gard this a reasonable choice for a relativitic ex- 
panding flow if the flow start from a Wien equi- 
lirium. This is seen in the behavior of the cosmic 
microwave background radiation in cosmology and 
the optically thick relativistic wind considered by 
GW98; the photon distribution remains a Planck- 
ian even after the optical thickness becomes much 
below unity. Anyhow, to obtain more exact pho- 
ton spectrum, we need to solve the radiative trans- 
fer coupled with dynamics, a challenging problem 
which will be pursued in future. Then, the ba- 
sic equations are the conservation laws of energy, 
momentum and numbers. For a steady and spher- 
ically symmetric flow, conservation laws of energy 
and momentum are given by 







and 



^ r [r\p + P)T^] + d 4=0, 



dr 



(1) 



(2) 



respectively. Here, p and P are the total energy 
density and pressure of electron-positron pairs and 
photons. Both components are described by ideal 
fluids and we use the units of c — 1. The bulk flow 
velocity and its Lorentz factor are denoted by j3 
and T, respectively. Momentum conservation law 
is alternatively written by 



1 dT 
fdr~ 



1 dP 

p + P dr 



= 0. 



(3) 



The energy density and pressure of electron- 
positron pairs are given by 



and 



p c ± = 2ra c n c (7) 



P c ± = 2m c n c 



(4) 



(5) 



respectively. Here, n c is the number density of 
electrons (=that of positrons) and the total num- 
ber density of leptons are 2n c , 9 is the temper- 
ature normalized by the electron mass m c and 
( 7 ) = K 3 (l/9)/K 2 (l/9) - 9 is the average Lorentz 
factor of electron thermal velocity, where Ki is the 
i-th order modified Bessel function of the second 
kind. 

As was mentioned above, we assume that pairs 
and photons are in a Wien equilibrium in the ini- 
tial state. In a Wien equilibrium, the energy den- 
sity and pressure of photons are given by 



and 



p 7 = 3m e n 7 # 



P 7 = m c n 7 9, 



(6) 



(7) 



respectively, where n 7 is the number density of 
photons. In a Wien equilibrium, n c and n 7 are 
related by 



n e K 2 {\/9) 
--^9^ = m 



(8) 



In numerical calculations of the flow we do not 
use eq.(8), but solve the following number conser- 
vation equations, assuming that pairs and photons 
take the same temperature, 



and 



(9) 



(10) 



where n e is the pair creation rate minus annihila- 
tion rate. The pair annihilation rate is approxi- 
mated by 



(«e)a 



= (Cannon;; 



= -a T n c 



(11) 



29 2 



Ln(1.120 + 1.3) 
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where v denotes the thermal velocity of pairs 
(Svensson 1982). The pair creation rate is given 

by 



fe) 2 


\K 2 {i/ey\ 




29 2 



(Ae)ann>(12) 



when the distribution function of photons is 
a Wien one (Svensson 1984). Since we ne- 
glect other photon production processes such 
as bremsstrahlung, there is a simple relation 

2h c = — hy. 

Then, we can solve for the bulk velocity, tem- 
perature and electron number density from the en- 
ergy, momentum and number conservation equa- 
tions given above. As the boundary conditions, we 
give the bulk velocity, temperature, and electron 
number density at the inner radius at r = r = 
2r g . The optical thick approximation is valid be- 
tween the inner radius and the photosphere. The 
radial coordinate of the photosphere r p h is defined 

by 

/?) = !. (13) 



/>OG 

/ dr 2n e cr T r(l 

Jr ph 



The optical thickness depends on the direction of 
the photon and for photons moving oblique to the 
radial direction the optical thickness is larger than 
that given above. But within the cone of 1/r 
around the radial direction the difference is small. 
Outside the cone, the optical thickness becomes 
much larger. Since we are treating relativistic out- 
flow, the most relevant is the radial optical thick- 
ness. 

In the optically thin region we treat pairs and 
photons separately. For pairs, energy, momentum 
and number conservation laws are given by 

^[r- 2 (Pe ± +Pe ± )r 2 /?]=F°, (14) 

^l[r 2 (Pe± + Pe±rP 2 ]+ ( ^ k =F 1 (15) 



and 



l__d 

r 2 dr 



[r 2 n c r/?] = n c 



(16) 



The right-hand side of energy and momentum 
conservation laws denotes the so-called radiative 
force, which is given through the Lorentz trans- 
formation of the radiative force in the fluid frame. 



The latter for Compton cooling is given by Phin- 
ney (1982) and Sikora et al. (1996) 



(Orad = -^(2n> T (( 7 2 - 1) - e,)T™ (] 7) 



and 



(F' 1 ) rad = (2n c )a T |i + ^( 7 2 )|7; 01 , (18) 

where e 7 denotes the phcnomenological Compton 
heating term described below and T^ v is the en- 
ergy momentum tensor of radiation field. The 
prime denotes the quantities in the fluid frame 
whereas those of the central black hole frame are 
represented without prime. Lorentz transforma- 
tion between these frames is described as 



F» = A^F'" and T^ v = A^A^T' 



/No- 



where 



A^„ 



r r/3 
173 r 



(19) 



(20) 



The Compton heating term is calculated as fol- 
lows. First, the average energy of photons in the 
fluid frame e is calculated and the equivalent pho- 
ton temperature # 7 = e/3 is defined. Noting that 
( 7 2 - 1) = 30K s (l/6)/K 2 (l/6), we put 



(21) 



The radiation field is assumed to be given by the 
free streaming approximation in the optically thin 
regime. The photons are emitted from the photo- 
sphere and they are relativistically beamed with 
the Lorentz factor of the photosphere. Although 
we take account of effects of photons on the dy- 
namics of electrons, we neglect the correspond- 
ing change of the photon field for simplicity. The 
equilibrium Lorentz factor of the radiation field is 
given by (Sikora et al. 1996) 



where 



r - = 7f( 1 + 7r=P 



1/2 



tX)0 _i_ fll ' 

7 ~ 7 



(22) 



(23) 



For F'°, we also include annihilation cooling 
term as 

(F' U n = 2m c ri e ( 7 > (24) 
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in addition to the Compton cooling and heat- 
ing term. We neglect other terms such as 
bremsstrahlung cooling, pair production from the 
photon field. Thus, our treatment is necessarily 
simplified to obtain detailed estimate of the final 
states of the flow, but we believe that the basic 
features of the flow are not much lost by these 
simplifications. 

3. Analytical solution 

In this section we describe analytic approxima- 
tion, which help to understand the physical situa- 
tion considered in this paper. 

3.1. Dynamics of Wien equilibrium pair 
plasma 

In the optically thick regime, adopting one fluid 
approximation and assuming that Wien equilib- 
rium (eq.8) is maintained in the flow too, we have 
two integrals from energy and number conser- 
vation laws. The energy conservation is writ- 
ten in the one fluid approximation from eqs. 
(1),(4),(5),(6),(7)&(8) 

Airr 2 m c n 7 {AO + 2f{6){6 + (7))} T 2 [3 = E, 

(25) 

where E is the total luminosity. Since the sum of 
the number of photons and pairs is assumed to be 
conserved, we obtain from eqs. (9)&(10) 



Airr 2 n 7 {l + 2f(0)}T0 = N, 



(26) 



where N is the total number flux of photons and 
leptons. 

Dividing eq.(25) by eq.(26), we obtain 



r = 



E 



1 + 2f(6) 



m e N A9 + 2f(8){8 + (7))' 



(27) 



In the limit of relativistic temperature, f(8) — > f 
and (7) — > 38, then we have 



in, 



E 1 
~NA8' 



(28) 



In the limit of non-relativistic temperature, f(6) 
and (7) — > 1, we have the same expression 



E 1 



m, 



,NA8' 



(29) 



It is to be noted that this expression is not an exact 
one at intermediate temperatures, in contrast to 
the case of complete thermal equilibrium(GW98). 
However, we use this as an approximate analytic 
solution in this section. This relation implies that 
the bulk Lorentz factor increases as the tempera- 
ture decreases, i.e., that the plasma is accelerated 
thermally in compensation of the internal energy. 
The factor E/m c N represents the mean energy 
of particles except for a numerical factor of order 
unity. Supposing that T at the boundary is order 
of unity, the bulk Lorentz factor attains an order 
of 10, if the plasma remains optically thick during 
an order of magnitude drop in the temperature. 
The temperature at the boundary needs not nec- 
essarily be highly relativistic as long as it is high 
enough to produce pairs. But, it is important to 
determine the ratio of pairs to photons at various 
radii. The actual kinetic power of pairs at infin- 
ity depends on the amount of surviving pairs too, 
which is considered in the following subsections. 

The photon number density is also represented 
by the temperature as (for T » f ) 



N 2 m c 40 + 2/(0)(0 + < 7 )) 

f) 

E 47rr 2 {1 + 2f(8)} 2 
At the relativistic temperatures, it becomes 



N 2 m A8 



7 E Airr 2 3 ' 
while at non-relativistic temperature, it is 

N 2 m 



E Airr 2 



A8. 



(30) 



(31) 



(32) 



The difference of the factor 3 means that at non- 
relativistic temperatures all the pairs annihilate 
and are converted to photons. Anyway, the total 
number flux of photons and electron pairs is con- 
served and total number density is approximated 
by 

N 2 m 

n 7 + 2n c = — C -A8. (33) 
7 E A-nr 2 

Under the approximations of T = E/Am e 8N 
and the above expression for n 7 +2n e , the equation 
of motion (eq. 3) becomes 



+ 



(l + 2f(8))r 2 



8dr 8[A9 + 2f(9){8+( 1 ))}dr y r 2 



d (S) = 0- (34) 
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Thus, both for relativistic and non-relativistic 
temperatures, we obtain 



9 oc r 



(35) 



Then, we obtain T oc r and n 7 + 2n c oc r~ 3 . Thus, 
thermal acceleration works for this Wien equilib- 
rium pair outflows. 

In order that optical thick condition is satis- 
fied, the optical thickness at the boundary should 
be at least larger than unity. We suppose that the 
boundary temperature 9 is at least mildly rela- 
tivistic and the number densities of photons and 
electrons are comparable (f(6o) — !)• Thus, we 
have 



n -f,0 — n cfi 



1 

3^7,1 

1 m e N\ 

1 E 
487rm c r%0 o T Q 2 



(n 7i0 + 2n C! o) 



(36) 



From eq.(13), the optical thickness at the bound- 
ary is estimated as 



TO 



n c n<TTro 



1 r g m p 1 



E 



(37) 



3r 72 r m c r 3 6» ^Edd 
where the Eddington luminosity is defined by 

L Edd ee 2^^ (38) 

and r g is the gravitational radius. Thus, in order 
to realize r 3> 1, the total luminosity should be 
at least an order of the Eddington luminosity for 
reasonable ranges of T and 9 at the boundary 
r = 2r g . 

The optical thick solution obtained in this sub- 
section is summarized as 



r = r - 

ro 
r 



and 



n 7 + 2n e = 3n c . ( — ) 
r 



(39) 
(40) 
(41) 



3.2. Photosphere and freeze-out radius 

The optical thickness decreases with radius and 
the flow eventually becomes optically thin. For 
relativistic temperatures, the optical thickness is 
given by 



n c a T r n c . a T r 



3r 



3r 



TO 



'•{) 



(42) 

Defining the radius of the photosphere r p h as the 
radius where r = 1, we obtain 



1/3 

r P h ^ V r . 
The temperature at the photosphere 6* p h is 



7ph 



1/3 



(43) 



(44) 



If 6> p h is smaller than unity, the real photosphcric 
radius is smaller than r p h given above because n e 
becomes smaller than n 1 . 

On the other hand, equating annihilation time 
scale with dynamical time scale, we obtain the 
freeze-out radius rf r as 



(n e ) 8 



r 



o T 



1/3 



ro 



1/3 



r ph . (45) 



These two radii are almost the same. Although 
the above expression for r rr is correct for 9 p h > 1, 
the coincidence of the two radii itself is valid even 
if 9 p h < 1. Thus, copious pairs exist in the opti- 
cally thick regime because of the photon photon 
pair production counterbalancing the pair anni- 
hilation, unless the temperature becomes too low 
at r p h- In the optically thin regime, a large frac- 
tion of pairs can survive because the annihilation 
time scale becomes longer than the dynamical time 
scale. 

The behavior of the freeze-out depends on the 
temperature at the photosphere 9 p h- For # p h < 1, 
the number density of pairs becomes small com- 
pared to the photon density but pairs can still 
carry a significant amount of kinetic power in the 
rest mass energy. For # p h > 1, the number density 
of pairs is comparable to the photon density and 
the kinetic power of pairs accounts for more than 
half of the total power. Surviving kinetic power, 
however, is affected by interactions with radiation 
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Fig. 1. — Photospheric temperature in the E — 9 a 
plane. The upper-left region of # p h = 1 cor- 
responds to the non-rclativistic freeze-out with 
p h < 1, while the lower-right region does to the 
relativistic freeze-out with 8 p h > 1. The lines of 
to =const. are also depicted. 
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Fig. 2. — Relation between 9 ann and # p h- The 
solution to equation (63) is shown. The dotted 
line represents 9 ann = 9 p h- 
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and 



(56) 



Since the electron number density follows the 
equation 



"e = -(o- a nnV) [n c 2 - n c . 

we obtain 

5n R ; 



hi ^ • 

(57) 



1 d 



2n ,cq (cr ann w) r 2 dr 



[n c , cq r 2 r/3] . (58) 



Since r oc r as long as the equilibrium is main- 
tained, we obtain 



Sn e w - 



r 



2r(cr a 



din rip 



rflnr 



(59) 



to the leading order. From equation (33), equi- 
librium number density of electrons is estimated 
as 



9m„N 2 



f(0) 



nr 2 E 2/(0) 
N 2 



rn, 



^exp(-i). (60) 



Then, equation (59) can be approximated as 



Sn e ~ - 



r 



2r(cr arin i>) 



2r(a ann v) 



n (61) 



for 9 « 1. 

We regard that the freeze-out occurs when 
5n c = n 0!eq is satisfied. We define the annihila- 
tion radius r ann as the radius where Sn c — n c . cq is 
satisfied to discriminate it from r& and r p h. The 
temperature at r an n is denoted by 9 ann . The equa- 
tion to determine r ann and 9 ann is given by 



,N 2 
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Recalling that r ann ~ r (r ann /r ) and r ann = 
To{9o/9 aml ) and using equations (36) and (37), we 
obtain the equation to determine 9 ann as 
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In Figure 2 the relation between # ann and 6> p h is 
shown and we see that the former is slightly higher 



than the latter for # p h < 0.6. The kinetic power of 
pairs estimated at the annihilation radius is given 

by 
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This simple estimate gives a somewhat higher 
value than that by equation (54). Thus, the ki- 
netic power of pairs is shown to be still fairly high 
and is comparable to the radiative power even for 
non-relativistic freeze-out. As for the Lorentz fac- 
tor at the photosphere, we expect that 



T p h 
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ph 
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The terminal Lorentz factor may be higher by a 
factor of a few owing to the additional thermal and 
radiative acceleration in the optically thin regime. 

4. Numerical Results 

In this section, we present the numerical re- 
sults which are compared with analytic predictions 
made in the previous section. Numerical calcula- 
tions are performed for a wide range of values of 
the total luminosity E and the boundary temper- 
ature #o at r = r = 2r g . We fix r = \/3/2; 
other boundary values t , n 7 ^ and n c ^ and the 
total number flux N are easily obtained through 
several relations described in the previous section. 
We solve eqs. (1),(3),(9)&(10) for optically thick 
regime and eqs. (14),(15)&(16) for optically thin 
regime. The radiation field in the optically thin 
regime is calculated in the free streaming approx- 
imation as was described in section 2. The photo- 
spheric radius which marks the boundary between 
these two regimes is determined iteratively. 

As a typical example of non-relativistic freeze- 
out, we show in Figure 3 the numerical results for 
E/L Edd = 24.9 and 9 Q = 2 with r = 82.4. Fig- 
ure 3(a) shows the behavior of the bulk velocity 
Tfi and temperature 9. It is seen that the pho- 
tosphere is located at r p h = 7.45r g and that the 
temperature decreases rapidly with radius and be- 
comes # p h = 0.41 at the photosphere. The bulk 
Lorentz factor at the photosphere turns out to be 
r p h/3 P h = 5.06. These numerical values reasonably 
agree with the analytic predictions of r p h = 8.7r ff , 
# P h = 0.46 and r ph = 5.3. In particular, it is 
seen that a factor of 5 decrease in temperature 
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matches very well a factor of 5 increase of the bulk 
Lorentz factor, which confirms the basic feature of 
the thermal acceleration of the pair plasma in the 
optically thick regime. The differences between 
analytic and numerical results may be due to some 
simplifications made in the analytic estimates such 
as r 3> 1 and pair annihilation process. 

It is particularly important that pair outflow is 
further accelerated outside the photosphere. The 
numerical result shows that the terminal Lorentz 
factor attains = 15.9, more than three times 
that at the photosphere. Most of this factor is ac- 
counted for by the internal energy of pairs at the 
photosphere with an additional effect of radiative 
acceleration. It should be noted that since even 
at 6 p h =0.41, 8 + (7) =2.3, pure thermal accel- 
eration can attain <~ 12. Moreover, beamed 
photons from the photosphere heat the pairs by 
Compton interaction since the effective tempera- 
ture of photons is larger than electron temperature 
outside the photosphere as is seen in this figure. 
So the thermal acceleration works efficiently while 
radiative acceleration also makes a small contribu- 
tion to bulk acceleration at larger distances. This 
additional bulk acceleration by a factor of 3 is 
an important element in comparison with obser- 
vations. 

Figure 3(b) shows the relative number fraction 
of pairs and photons as a function of the radial 
coordinate. As is seen, the pair fraction decreases 
steadily according to the decrease of temperature 
and annihilation continues to 2r p h- Thus, the an- 
alytic estimate of rf r and r ann may be uncertain 
within a factor of 2. But, the degree of the de- 
crease is modest and about 20% of pairs annihilate 
outside the photosphere. To sum up, about 60% of 
the initial pairs survive in this case. This confirms 
our prediction that the pair annihilation problem 
can be avoided when we consider an initial Wicn 
equilibrium state. 

Figure 3(c) shows the luminosity of pairs and 
radiation as a function of the radial coordinate. 
The total luminosity turns out to remain constant 
function of radius as it should be. Inside 
the photosphere, the kinetic power of pairs slightly 
decreases with radius, which is because pairs are 
converted to photons as the temperature decreases 
keeping a Wien equilibrium. The kinetic power of 
pairs is kept almost constant outside the photo- 
sphere as it should be for bulk acceleration by their 
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Fig. 3. — (a) Behavior of velocity and tempera- 
ture for E/L E dd = 24.9 , 9 = 2 with r = 82.4. 
The solid line denotes the bulk Lorentz factor of 
the pairs T(3, and the dash-dotted line denotes 
the equilibrium Lorentz factor of radiation field 
r eq /3eq- The dashed line denotes the temperature 
of pairs, and dotted line denotes the equivalent 
temperature of radiation outside the photosphere, 
(b) Fraction of particle and photon numbers for 
E/L Edd = 24.9, O = 2 with t = 82.4. The solid 
and dashed lines denote the number fractions of 
pairs and photons, respectively, (c) Luminosities 
of various components for E/L^dd — 24.9 , 9 = 2 
with ro = 82.4. The dashed line denotes the total 
luminosity of pairs and radiation. The solid line 
denotes the kinetic power of pairs while the dot- 
ted line denotes the power carried in a form of the 
rest mass of pairs. Thus, the interval between the 
dashed and solid lines denotes the luminosity of 
radiation and that between the solid and dotted 
lines denotes the power carried in a form of the 
thermal energy of pairs. 



own thermal energy. A small change is due to the 
interaction with photons. But, it is seen that the 
energy gained from Compton interaction almost 
cancels out with a small but finite amount of an- 
nihilation cooling. The terminal kinetic power of 
pairs turns out to be 61% of the total luminosity, 
while remaining 39% is carried away by photons 
emitted from the photosphere. Since the matter 
at the photosphere is moving at T = 5, distant 
observers see the radiation peaked at an energy 
of r# p h ~ 2, i.e., around lMeV. Thus, this model 
predicts that strong MeV emission should accom- 
pany the relativistic jet formation. 

We have also calculated for various boundary 
values and the results are tabulated in Table 1. 
As is seen, the terminal Lorentz factor proves to 
be 8 ~ 60 and the terminal kinetic power of pairs 
accounts for 50 ~ 80% of the total luminosity for 
1 < O < 5 and 0.8 < E/L Edd < 80. Thus, the 
Wicn fireball turns out to be quite an efficient 
mechanism for the production and bulk acceler- 
ation of relativistic jets. As for the dependence of 
the physical quantities at the photosphere, they 
reasonably agree with the analytic predictions of 

-V3 a , „ p-l/304/3 



oc r 1/3 oc EV*9~ 



' ph oc E-V 6 6l'° and 



r p h oc E ' O ■ Figure 4 shows r p h as a func- 
tion of To and figure 5 does p h as a function of 
t and then, we see that the analytic prediction 
agrees with numerical results, although the ana- 
lytic prediction is slightly larger than the numeri- 
cal result. 

Since the enthalpy of pairs at the photosphere is 
large, pairs are thermally accelerated outside the 
photosphere. The terminal Lorentz factor is thus 
basically determined by the boundary tempera- 
ture for relativistic freeze-out. For non- relativistic 
freeze-out, additional effect by radiation works to 
further increase the terminal Lorentz factor. This 
behavior is clearly seen in Figure 6. The termi- 
nal Lorentz factor is above the naive estimate 40 
to a varying degree according to a varying degree 
of importance of additional thermal and radiative 
acceleration. It is to be noted that fairly large 
Lorentz factor can be obtained even for the initial 
temperature as low as 0.5. 

Finally, in Figure 7 is shown the efficiency of the 
outflow production, i.e., the ratio of the terminal 
kinetic power of pairs to the total luminosity. It is 
basically determined by the photospheric temper- 
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Fig. 4. — Relation between photospheric radius 
r p h and the initial optical thickness To- The solid 
line represents the analytic estimate r p h/Vo = 
To 1 / 3 . Open circles, solid diamonds, rectangles, 
triangles, crosses and stars refer to E/L Edd = 
83.1,24.9,8.31,2.49,0.831 and 0.249, respectively. 
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Fig. 5. — Relation between photo-spheric tem- 
perature # p h and the initial optical thickness to. 
The solid line represents the analytic estimate 
Oph/Oo = tq^ 1 / 3 . Open circles, solid diamonds, 
rectangles, triangles, crosses and stars refer to 
E/L Edd = 83.1,24.9,8.31,2.49,0.831 and 0.249, 
respectively. 



Fig. 6. — Terminal Lorentz factor Too for various 
boundary values. The abscissa is the temperature 
at the boundary 9q- The symbols are the same as 
Fig. 4. The solid line represents a naive expecta- 
tion of Too = 4# - 
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ature. For relativistic freeze-out, the efficiency is 
about 10% larger than the canonical value of 2/3. 
This excess is due to additional radiative accelera- 
tion. For non- relativistic freeze-out, the efficiency 
becomes small because pair density decreases ow- 
ing to pair annihilation. At 8 p h = 0.2 the efficiency 
is about 30% and it will be very small for lower 

Summarizing these results, we have found that 
electron-positron pairs can be thermally acceler- 
ated up to the bulk Lorentz factor of more than 10 
and that their kinetic power becomes comparable 
to the total power, provided that Wien equilibrium 
states of pure electron-positron pairs are prepared 
at the inner boundary with a temperature of a few 
times electron mass and the optical thickness to 
the scattering of more than around 5. The corre- 
sponding total luminosity exceeds the Eddington 
luminosity by a factor of more than unity in the 
spherical symmetry. If we assume that the outflow 
is conical and that the opening angle of the outflow 
covers f2 steradian, the total power is smaller by 
a factor of Air/fl. Since is typically as small as 
10~ 2 , the required total luminosity can be as small 
as 1% of the Eddington luminosity. The mecha- 
nisms of the realization of pure electron-positron 
pairs and initial collimation of the outflow remain 
to be a theoretical challenge. However, we believe 
that this is the most successful specific mechanism 
of the bulk acceleration of relativistic outflow up 
to the Lorentz factor of 10 with high efficiency. 
This Wien fireball model predicts that radiation 
from the photosphere should appear as a strong 
emission around a few MeV with a comparable lu- 
minosity to the jet kinetic power. 

5. Discussion 

Finally, we discuss several issues related to the 
Wien fireball model. The first problem is the re- 
alization of the initial state. Most probably such 
a state is prepared in relation to the hot accre- 
tion disks or hot accretion disk corona. Previous 
studies on pair concentration in accretion disks 
have shown that the pair concentration in the disk 
is not so large under the steady state condition 
and that the steady state solution does not exist 
when the luminosity is higher than about a few 
percent of the Eddington luminosity (Kusunose 
and Takahara 1988). This is due to the ineffi- 



Fig.7 




Fig. 7. — Efficiency of the outflow production; i.e., 
the ratio of the terminal kinetic power of pairs to 
the total luminosity. The abscissa is the photo- 
spheric temperature. The symbols are the same 
as in Fig. 5. The dotted line refers to the effi- 
ciency of 2/3. 
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ciency of the pair annihilation and suggests that 
pair runaway occurs for such high luminosity. In- 
stead of the pair equilibrium condition, if we al- 
low free escape of pairs from the disk, Yamasaki, 
Takahara, and Kusunose (1999) have shown that 
much higher pair concentration is obtained and 
that outflow luminosity of pairs can be very large 
based on a simple analysis. This seems to be a 
promising way to produce the initial state of the 
present outflow model, since Yamasaki, Takahara, 
and Kusunose (1999) predicts the optical thick- 
ness of pairs is around a few and the temperature 
of a few. Although pure pair escape was assumed 
in Yamasaki, Takahara, and Kusunose (1999), in 
reality baryon load should be carefully examined. 
If the baryon density is more than 1% of the pair 
density, the terminal Lorentz factor will be signifi- 
cantly decreased. This is a common problem with 
the thermal fireball for GRBs. 

As for the photon field, we assumed that Wien 
equilibrium photons are prepared at the bound- 
ary and subsequent photon production process is 
neglected except for pair annihilation. If the ini- 
tial optical thickness is so large, this seems to be 
a fair assumption but for optically thickness of 
a few to 10, this is not true and photon spec- 
trum will not be a Wien one but near the Comp- 
tonized bremsstrahlung. Moreover, in order to 
have a large pair concentration, the amount of 
soft photons must be small because they will de- 
crease the electron temperature and pair produc- 
tion. In any case, radiative transfer including such 
as bremsstrahlung and Compton scattering should 
be examined in a future work. This is also im- 
portant to predict the accompanying photon spec- 
trum with jet production. It is to be noted that 
whole process should be treated dynamically, since 
the pair escape is the most critical assumption for 
this mechanism. 

Our calculation indicates that super Eddington 
luminosity is needed for an relativistic jet produc- 
tion provided that the flow is spherically symmet- 
ric. As was mentioned in the previous section, 
if some collimation mechanisms exist near the 
boundary, the required luminosity will be smaller 
and the real luminosity can be sub-Eddington. 

For many years, radiative bulk acceleration 
mechanism has been regarded ineffective in pro- 
ducing relativistic jets. But most of them treat 
optically thin case (Phinney 1982; Inoue and Taka- 



hara 1996). Since the equilibrium Lorentz factor 
of photons emitted from typical accretion disks is 
not large enough near the accretion disk, radiative 
force acts as radiation drag rather than as accel- 
eration. So the attainable Lorentz factor of bulk 
flow is at most about 3. This is because radia- 
tion field is not collimatcd enough to the jet di- 
rection. Further, for optically thin plasma, only 
a small fraction of the radiation power is con- 
verted to plasma kinetic power. Thus, the op- 
tically thin radiative acceleration is ineffective in 
general. In contrast, our present model assumes 
that the pair plasma starts to flow out with op- 
tically thick condition. This means that photons 
move at the same speed with pairs since they are 
strongly coupled and that they can be thermally 
accelerated. Although the pair outflow becomes 
optically thin outside the photosphere, radiative 
acceleration acts effectively because radiation has 
been collimated at the photosphere and the equi- 
librium Lorentz factor becomes large enough. Of 
course, the radiation drag from ambient photons 
will somewhat decrease the bulk velocity of the 
jets, but the effect will be small since the beamed 
radiation field dominates over ambient photons for 
typical cases. 

Another argument against pair jets has been 
the annihilation problem (Celotti and Fabian 
1993; Blandford and Levinson 1995). When a 
large amount of cold pairs exist in a compact re- 
gion, they will inevitably annihilate before they 
escape. Since in our model pairs are hot and ac- 
companied with high energy photons, pair produc- 
tion works as well. Inside the photosphere, pair 
annihilation is balanced with pair production and 
outside the photosphere annihilation time scale 
becomes longer than the dynamical time scale. 
Thus, the most pairs can survive to infinity and 
can avoid the annihilation problem with pair jets. 

This work is supported in part by the Grant- 
in-Aid for Scientific Research of the Ministry of 
Education and Science No. 11640236 and by Re- 
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Table 1: Numerical results for various boundary 
conditions. First, second and third columns refer 
to the total luminosity, boundary temperature and 
boundary optical thickness. The 4-th through 6-th 
columns show the numerical results for the quan- 
tities at the photosphere, while 7-th through 9-th 
columns show the numerical results for those at 
infinity. 
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